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Abstract Aerosol backscatter coefficient data were examined from two flights near Japan and 
Hawaii undertaken during NASA’s Global Backscatter Experiment (GLOBE) in May June 1990. 
During each of these two flights the aircraft traversed different altitudes within a region of the 
atmosphere defined by the same set of latitude and longitude coordinates. This provided an ideal 
opportunity to allow flight level focused continuous wave (CW) lidar backscatter measured at 
9. 1 1-pxn wavelength and modeled aerosol backscatter from two aerosol optical counters to be 
compared with pulsed lidar aerosol backscatter data at 1 .06- and 9.25-prn wavelengths. The best 
agreement between all sensors was found in the altitude region below 7 km, where backscatter 
values were moderately high at all three wavelengths. Above this altitude the pulsed lidar 
backscatter data at 1 .06- and 9.25-pm wavelengths were higher than the flight level data obtained 
from the CW lidar or derived from the optical counters, suggesting sample volume effects were 
responsible for this. Aerosol microphysics analysis of data near Japan revealed a strong sea-salt 
aerosol plume extending upward from the marine boundary layer. On the basis of sample volume 
differences, it was found that large particles were of different composition compared with the 
small particles for low backscatter conditions. 


1. Introduction 

The study of tropospheric aerosols in remote regions has 
now become an important component for modeling cloud 
physics, radiative transfer, global transport, and climate forc- 
ing effects in the atmosphere [Coakley and Cess , 1985; Al- 
brecht , 1989; Charlson et al. t 1992], Often these investiga- 
tions are done in the remote troposphere using moving plat- 
forms such as aircraft with sensors, which collect data on aero- 
sol properties at aircraft flight level and at a distance from the 
aircraft. Sensors used in aerosol studies vary widely in the way 
the aerosol characteristics are obtained. In the case of aerosol 
backscatter coefficient measurements these data can be (1) 
measured directly with range- resolved pulsed lidars focused at 
infinity, (2) derived from aerosol microphysics at flight level, 
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or (3) measured directly at aircraft level with a focused continu- 
ous wave (CW) lidar. Typically, pulsed lidars are unable to 
measure aerosol backscatter within - I km of the aircraft (due to 
nonoverlap of transmitted and received beams as well as op- 
tomechanical effects). This precludes directly comparing flight 
level derived aerosol backscatter with data measured by these 
lidars. Therefore one cannot compare the effects of lidar pulse 
length (i.e., sample volume) and time averaging with flight 
level backscatter data viewed horizontally from the same spa- 
tial region at the same time. The effect of sample volume is 
particularly important because of the large sampling volumes 
associated with pulsed lidars compared with sensors used to 
measure aerosol size distributions at flight level, i.e., lidars 
may more thoroughly sample sparsely distributed large parti- 
cles. 

The National Aeronautics and Space Administration (NASA) 
conducted the Global Backscatter Experiment (GLOBE II) in 
May-June 1990 [Bowdle et al., 19911 over and around the re- 
mote Pacific Ocean. The centerpiece of GLOBE was an instru- 
mented NASA DC-8 aircraft that operated in the troposphere. 
The aircraft carried a suite of sensors which allowed aerosol 
backscatter coefficients to be derived remotely using pulse li- 
dars and at flight level using CW lidar. Several of the flights 
were arranged such that the aircraft traversed the same region of 
atmosphere which had shortly beforehand been remotely sam- 
pled by two pulsed lidar systems. In this paper the results of 
aerosol backscatter coefficient intercomparisons at several 
wavelengths are presented, using data from five different sen- 
sors. Explanations are provided for aerosol backscatter coeffi- 
cient variations in terms of instrument performance, aerosol 
microphysics, and composition prevailing at the time. 
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2. Flight Details 

Aircraft flight details and meteorological conditions are 
presented for two flights which took place near Japan and Ha- 
waii as part of the GLOBE 11 Pacific Ocean survey mission us- 
ing the NASA DC-8 aircraft. 

2,1. Japan Flight 

This flight (number 13 in the GLOBE II flight scries) look 
place on June 3, 1990, from 0108 to 0808 UTC, over the Pa- 
cific Ocean, SE of Tokyo, Japan. The data reported here were 
obtained from 0240 to 0545 UTC for all sensors. This portion 
of the flight basically consisted of a series of -280-km level 
flight legs from NW to SE at 1 1.5, 9.8, 6.5, and 3.2 km be- 
tween (26.10° N, 146.70° E) and (23.74° N, 148.66° E). Each 
level flight leg (sec Figure la) was followed by a gradual de- 
scent in the opposite direction to the starting point of the next 
level flight leg. The slow descents provided detailed vertical 
profiles from the flight level instruments, for comparison with 
the pulsed lidar profiles obtained remotely along the same 
flight track. An exception to this pattern was the rapid descent 
from 1 1.5 km where the DC-8 reached the 9.8-km level well 
before the NW end point of the intercomparison ground track. 
In this case the aircraft retraced a large portion of this flight 
level twice but in opposite directions. 

Meteorological conditions during this flight featured a 40 
m s 1 westerly jet above 8-km altitude near Tokyo. The jet 
gradually weakened as the aircraft progressed southeastward 
during the first 1 .5 hours of the flight. Considerable cloudi- 
ness was found throughout the troposphere, especially near 



Figure 1. Flight tracks near (a) Japan on June 3, I°90, and (b) 
big island of Hawaii on May 18-19, 1990. The flight southeast 
of Tokyo consisted of four stacked flight levels between 
points A and B at the altitudes indicated in the insert, while the 
Hawaii flight was composed of four closed loops (at constant 
altitudes indicated). 


Japan, associated with a cold front trailing from a midlatitude 
cyclone northeast of Japan. Additional cyclogenesis was evi- 
dent southwest of Kyushu in response to the approach of a new 
synoptic wave moving out of Mongolia. Between 0230 and 
0300 UTC, as the aircraft approached its farthest point from 
Tokyo, it encountered weak easterlies near 11-km altitude, in- 
dicating the aircraft's traversal of the core of the belt of sub- 
tropical anticyclones south of Japan. The aircraft then exe- 
cuted its first U-turn and began its series of flight legs at ever 
decreasing altitudes. During these descending flight legs, 
winds at flight level shifted to the northwest quadrant, al 
speeds of 10 m s' 1 or less (see Figure 2) and remained from that 
quadrant until the aircraft ultimately regained altitude and 
headed back toward Tokyo. At no lime did the aircraft encoun- 
ter clouds during the maneuver. However, a complex pattern of 
haze layers and sudden changes in ambient dewpoint at some 
flight levels during the descending flight legs were observed. 
Examining Figure 2a shows that the wind direction during the 
highest stepped descents (9.8 and 11.5 km) was mostly at an 
angle of 50" to the flight track (these angles arc indicated by 
the two dashed lines in the figure). With wind speeds between 5 
and 10ms 1 the aircraft exhaust would be sufficiently removed 
from the area when the lidar remotely sensed the higher flight 
levels from below. Contamination at the higher altitudes 
would be of concern because of low aerosol backscatter gener- 
ally prevailing. 

2.2. Hawaii Flight 

This flight (number 5 in the GLOBE II flight scries) look 
place around the big island of Hawaii on May 18-19, 1990, 
from 21 14 to 0142 UTC. Figure lb shows a three-dimensional 
(3-D) view of the flight track above the island of Hawaii. It 
consisted, in part, of four complete loops near the coast of the 
island at 4.0-, 6.4-, 9.4-, and 12.5-km heights, maintaining 
the same course position for each loop and rapidly climbing to 
the next level at the NW corner of the island. 

Meteorological conditions at the time of this flight were in- 
fluenced by the recent passage of a sharp middle and upper tro- 
pospheric trough, the trailing portions of which were just 
northeast of the islands, and by a convective cluster along the 
Jntertropical Convergence Zone (ITCZ) south of the islands. 
Lower tropospheric winds were generally from the northeast al 
10 m s' 1 backing to the north at roughly 20 m s' J in the middle 
troposphere. Near 12-km altitude, a westerly jet having speeds 
from 40-50 m s 1 prevailed. The atmosphere was moist below 
700 hPa, and surface dewpoints near 20 n C sustained scattered 
convective showers topped by layers of detrained altocumulus. 
A thermal inversion was present near 700 hPa, and substantial 
amounts of haze were trapped below this layer. Above 700 hPa, 
the air was much drier, with only patches of cirrus near the is- 
lands. A region of more widespread cirrus detrained from the 
ITCZ was visible to the south. 


3. Measurements 

The lidar instrumentation consisted of three lidar systems 
operating at six wavelengths. Pulsed lidar systems were pro- 
vided by NASA Jet Propulsion Laboratory (JPL) operating at a 
single wavelength of 9.25 Jim [Menzies and Tran , 1994] and 
by NASA Goddard Space Flight Center (GSFC) operating at 
0.53-, 1.06- and 1.54-pm wavelengths [Spin hi me et al., 
1997]. Two focused CW Doppler lidar systems operating at 
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9. 1 I - and 10.59-pm wavelengths were provided by NASA Mar- 
shall Space Flight Center (MSFC) [Rothermel et al., 1996]. All 
three lidars were calibrated to measure absolute aerosol 
backscatter coefficient ((3, trf 1 sr' J ) as described in the reference 
to each system listed above. 

The aerosol size distribution at flight level was measured us- 
ing two different optical particle counters (OPC). A wing- 
pod-mounted forward scattering spectrometer probe (FSSP), 
model 300, from NASA Ames Research Center (ARC) covered 
the diameter range 0.35-20 jam with 10-s temporal resolution 
[Pueschel et al. , 1994]. The University of Hawaii laser optical 
particle counter (LOPC) covered the size range 0.16- to 7.0- pm 
diameter using sampling times that were of the order of several 
minutes [ Clarke , 1991, 1993]. A detailed summary of these 
sensors and their sampling characteristics has been given 
elsewhere [Cutten et aL, 1996]. In addition, an ultrafine con- 
densation nuclei (UCN) counter measured nuclei between 3- and 
15-nm diameter, while a second counter measured condensation 
nuclei (CN) above 15-nm diameter f Clarke , 1993]. The latter 
counter also measured nuclei in ambient and heated aerosol to 
provide data on the CN refractory component. Table 1 lists the 
uncertainty in (3 measured by each of the lidars and the uncer- 
tainty in the OPC-measured aerosol size distributions. The fi- 
nal uncertainty in the OPC-modeled (3 will be affected by un- 
certainties in the modeled aerosol composition, which are dif- 
ficult to quantify. However, data indicate that the LOPC- 
modeied P uncertainty is of similar magnitude as the CW lidar 

p. 

3.1. Size Distribution Data 

The FSSP size distribution data were processed as described 
by Cutten et al. [1996]. This consisted initially of resizing the 
FSSP sample bins using LOPC-modeled aerosol refractive in- 
dex data at 0.55-pm wavelength. The \0-s cfn/d\ogD it data were 
then averaged over 100 s and fitted to single-mode lognormal 
curves. The fitting process minimizes the effects of poor count 
statistics in the large particle regime. The FSSP underestimates 
bin counts in the 0.8- to 3-pm diameter range because of inade- 
quate sampling; this size range dominates infrared P 
[Srivastava et al ., 1992]. Consequently, it has been found 
[Cutten et al., 1996] that FSSP-modeled P less than !0' 1(1 nV 1 sr* 

1 in the infrared wavelength region (9-1 1 pm) can be as much 


as a factor of 3 smaller than measured data even with the curve 
fitting applied. 

FSSP-measured size distributions were used to calculate P 
from Lorentz-Mie theory and LOPC-modeled aerosol refractive 
index data. The refractive index data were derived from inferred 
mixtures of sulfuric acid, ammonium sulfate, and dust using 
volumetric mixing rules. A growth correction was included due 
to relative humidity differences between ambient and LOPC 
sampling volume. Calculations for P were done at 1.06- and 
9.1 1-pm wavelengths. Since the differences in refractive index 
of acid sulfate at 9.11- and 9.25-pm wavelengths are very 
small, measured and modeled p at 9.1 1-pm wavelength can be 
directly compared with measured p at 9.25-pm wavelength. 
Here p was also modeled from LOPC size distribution data 
based on aerosol mixture models described by Srivastava et aL 
[1997] that were most suitable for different aerosol loading 
conditions. Using the differences between the thermally seg- 
regated LOPC size distribution, volumes of different composi- 
tions were inferred. In cases of midtroposphere clean aerosol 
loading conditions, the mixed phase sulfate model was used to 
calculate refractive indexes of each LOPC size bin. For higher 
loading cases with a dominant dust component, a composite 
model was used for estimating the refractive index of the aero- 
sols. These were used in Lorentz-Mie calculations to get P for 
different wavelengths. 

3.2. Lidar data 

JPL and GSFC pulsed lidar P data were processed in two 
ways. Data were extracted from those portions of nadir or ze- 
nith profiles within 0.5-km altitude of any selected flight 
level. These backscatter data were then averaged, either in alti- 
tude for time series along each leg or in time and altitude for a 
composite on each leg. These averaged data could then be 
compared directly with modeled or measured p at flight level, 
either in a time series format or in an altitude versus p format. 
The second way of processing the data involved temporally 
averaging the measured profiles over time during each con- 
stant-altitude leg from four flight levels for the Japan flight 
and from two flight levels for the Hawaii flight. These profiles 
were used to evaluate the range correction profile for each 
pulsed lidar. They also provided an indication of any temporal 
variability in the aerosol field during each intercomparison 
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Figure 2, Time series of (a) wind direction, (b) w'ind speed, and 
(c) relative humidity along the flight track near Japan (Figure 
la). The dashed lines in Figure 2a show track angles for the in- 
tercomparison period. Dashed line in Figure 2b is aircraft alti- 
tude. 


experiment. Lidar (3 profiles were not always available at each 
level in the two case studies. This was because the iidars were 
at times pointing in the zenith and not nadir direction when 
the aircraft was at the higher altitudes, problems with the in- 
strumentation, or backscatter from intervening clouds. Aver- 
aging processes excluded any GSFC or JPL P data contami- 
nated by cloud returns. 


In the case of the 9.11- and 10.59-pm P data from the 
MSFC CW Doppler lidar, all 3-s P data falling within the 
FSSP-modeled P 100-s averaging period were averaged to pro- 
vide coregistered data. For the time series plots, a 300-s 
“boxcar” running average was applied to both sets of data. 


4. Interconiparison of OPC Modeled 
and Lidar (3 

This section will examine and evaluate the iniercomparison 
of measured and modeled p data from pulsed Iidars and in situ 
sensors for flights near Japan (section 4.1) and Hawaii 
(section 4.2). 

4.1. Japan Flight 

Before beginning the p intercomparison discussion, the 
dominant aerosol feature observed in this flight will be pre- 
sented. Plate 1 shows a time-altitude false color cross section 
of logarithm of P 106 for the whole of Japan flight 
(subsequently referred to as 1.06 -jam colorized P plot). The 
near-continuous, high-resolution (60 s) vertical profiles mak- 
ing up this image clearly show an aerosol plume extending 
vertically from the marine boundary layer (MBL). This feature 
is clearly reproduced when the aircraft repeated its track at dif- 
ferent altitudes. The corresponding 9.25-pm colorized P image 
(not shown) did not show this feature as clearly because (1) 
most of the backscatter is from larger particles, and (2) a 
longer profile integration time (100 s) was used along with a 
20-s gap between each profile. 

4.1.1. L06 pm. A time series plot of measured and modeled 
p (()f> for the two lowest altitudes is shown in Figure 3. At both 
altitudes, the time series for each instrument shows transects 
of the plume denoted by regions A, B, and D in Plate 1. The 
plume provided a prominent intercomparison target, with P 
values at each altitude ranging over nearly an order of magni- 
tude along each -280-km transect. The highest values, -4x10 K 
m 1 sr'\ occur in the plume peak at 3.2-km altitude. The lowest 
values, ~7x 1 0 y m' 1 sr ‘, occur in the plume fringes at 6.5-km 
altitude. In all cases, P, ()6 values are well above the nominal de- 
tection threshold of the GSFC lidar and the count-limited detec- 
tion thresholds of the FSSP and LOPC for modeled P, w> [Cutten 
ei c d . , 1996]. 

At both altitudes, plume transects for each instrument show 
well-defined peak and fringe features. Details vary with alti- 
tude, instrument, and sample period. In particular, as the flight 
progressed, the peak value occurred later in the transect, -4.8 
min (66 km) at 6.5-km altitude from 0240 to 0532 UTC and 
-2.4 min (33 km) at 3.2-km altitude from 0240 to 0337 UTC. 
(DC-8 true air speed was 230 m s' 1 .) These displacements are 
consistent with plume advection by the mean along-track wind 
at each altitude (50 and 34 km at 6.5- and 3.2-km altitude, re- 
spectively). Subsequent discussions of the intercomparison re- 
sults at these two altitudes will assume that plume properties at 
each altitude did not change with the movement of the plume. 
This assumption allows comparisons of the magnitudes and 
temporal trends in the corresponding time series. 

The GSFC lidar data show excellent agreement (within 5%) 
in both magnitude and trend for almost all segments of the suc- 
cessive time series at both altitudes. The largest discrepancies 
occur in the last run, where the first few minutes show higher P 
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Table 1. Uncertainties for Lidar and Optical Particle Counter Data 


Instrument 

Wavelength, 

pm 

Sample 
Volume, nP 

Threshold (3 Level 
m' 1 sr' ( 

AP, 

m' 1 sr' 1 

An, 
cm -3 

GSFC pulsed lidar 

0.53 

4x10' 

5x10" 

11x10 s 

- 

GSFC pulsed lidar 

L06 

4xl0 5 

1x10" 

±1x10" 

- 

MSFC CW lidar 

9.11 

6xl0 2 

6x1 0 1 ’ 

+2.4x1 0 13 " 

- 

JPL pulsed lidar 

9.25 

6x10* 

1x10 11 

+2x10" 

- 




at 2 km from aircraft 



MSFC CW lidar 

10.59 

6xl0 3 

2xl0' n 

±7.2x10 13 * 


FSSP 

- 

1x10 3 

- 

- 

<±30% f 

LOPC 

- 

1.7x10’ 

- 

- 

<±20% 


* Note that the two CW lidar uncertainty values given in Table 2 by Cullen el til. [1996] were incorrect 
and should be replaced with these values. 

r FSSP uncertainty based on count statistics for size bins >0 8-pm diameter. 


by -50% and the last few minutes show lower P by -30% than 
the earlier runs. Both time series for the FSSP-modeied p [0(S 
match the trends in the corresponding GSFC time series 
(typically within 20%). However, at both altitudes the FSSP 
data are systematically lower than the GSFC data (a factor of 2 
at 3.2 km and a factor of 4 at 6.5 km), which is probably due to 
different sampling volumes of the sensors. By contrast, the 
peak and minimum values in the LOPC P, 06 time series show 
good agreement (within 10%) with the corresponding extrema 
in the GSFC data but significantly more disagreement (20- 
100%) in details of the temporal trends. In particular, LOPC 
peak values typically occur -5-7 min after FSSP peak values in 
the same time scries and therefore -5-7 min after the corre- 
sponding GSFC peak values. 

The plume properties at 3.2 km and 6.5 km were essentially 
invariant throughout the intercomparison period. They show 
interna] consistency in the GSFC p l()6 profiles for several in- 
dependent ranges and beam orientations (8.2-km nadir versus 
6.6-km nadir for the 3.2-km data and 5-km nadir" versus' 3.3-km 
nadir versus 3.3-km zenith for the 6.5-km data). They also 
confirm the inherent reliability of p J06 modeling using aerosol 
microphysical measurements from the LOPC, after standard 
corrections for sampling effects. Since the FSSP and GSFC 
trends showed good agreement, the discrepancies between the 
LOPC and GSFC trends can be attributed to Tong LOPC integra- 
tion times (-10 min, versus 100 s for the FSSP and 60 s for 
GSFC data) in the presence of strong systematic aerosol gradi- 
ents. 

The lime-averaged GSFC aerosol (3 profiles (Figure 4) re- 
veal very good consistency among the individual profiles con- 
sidering that they were measured several hours apart, with the 
biggest deviations occurring only at 4.7, 7, and 9 Tim There- 
fore aerosol temporal variability over the selected measure- 
ment periods can be considered to be minimal for the 1.06-pm 
wavelength throughout the troposphere. In fact, the residual 
discrepancies are due primarily to sampling successively dif- 
ferent portions of the same features as those features drifted 
across the intercomparison region. The overall trend in all of 
the profiles is to a lower (3 up to 10 km, which then remains 
constant at around I x I0 K m*' sr' 1 up to 15 km. Figure 4 also 


shows data for the GSFC lidar and FSSP and LOPC sensors in a 
different format at the four intercomparison altitudes. Here the 
time averaged P and corresponding upper and lower limits for 
each set of GSFC, FSSP, and LOPC data are plotted to indicate 
how the averaged data compare: Data used to compute the GSFC 
average p are the ±0.5-km data from all time periods available 
at that altitude. For the lower two altitudes the average and 
range for the GSFC, LOPC, and FSSP data all agree within a 
factor of 2-3. As noted previously, the GSFC and LOPC data 
show the best agreement, with the FSSP consistently showing 
lower values. At the two higher altitudes the agreement is gen- 
erally poorer, with the FSSP data again showing the largest 
difference. The LOPC data are in good agreement with the 
GSFC lidar data at 9.8-km altitude, while at the 1 L5-km alti- 
tude, the LOPC data are also lower. At these altitudes a pre- 
dominant fraction of the p in these cleaner regions was found 
to arise from a large number of smaller particles (<0.5-pm di- 
ameter) and an absence of larger particles. Hence, since the 
FSSP samples particles in a very restricted size range, the 
FSSP-modeled p will be underestimated. 


4.L2. 9,25 pm. The JPL lidar P data obtained at constant 
altitude during the Japan flight arc shown as a lime series in 
Figure 5, again at the same two flight altitudes as for the 1 .06- 
pm data (Figure 3). Along with these data are also shown the 
100-s FSSP-modeled, LOPC-modeled, and 100-s CW lidar 
measured p, all at flight level. Figure 6 shows time-averaged 
JPL profiles for each of the three time intervals examined, 
along with a set of averaged data from each sensor. At the two 
lower altitudes (3.2 and 6.5 km; Figure 6), mean P values agree 
to within a factor of -3 for all four flight level sensors. How- 
ever, the spread for the CW lidar data is quite large at the 6.5- 
km level, as well as for the LOPC data at 3.2 km. At 9.8-km al- 
titude (not shown in Figure 5) the mean values differ by over a 
decade. However, the CW P lidar data and FSSP- and LOPC- 
modeled p agree to within a factor of -1.5 at 1 1.5-km level. 
The first 3 km of the 0518-0538 UTC profile up to 7 km has a 
similar shape to the corresponding portion of the 0424-0442 
UTC profile although displaced to higher P values. The profile 
was truncated at 9 km because of excess noise between 9 and 
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Plate 1. False-color cross section of logarithm of 1.06-pm aerosol backscatter near Japan on June 3, 1990. 
Regions marked A, B, C, and D are periods used for sensor intercomparison. 
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Figure 3. Time series plot of measured and modeled 1,06-jim aerosol backscatter coefficient p for selected 
periods from the Japan flight on June 3, 1990. Lidar data are averaged over ±0.5 km of the two flight alti- 
tudes studied. The p data modeled from the forward scattering spectrometer probe (FSSP) and laser optical 
particle counter (LOPC) flight level sensors are superimposed. Dashed horizontal line at 1 x 10' y m' 1 sr’ 1 
represents the Goddard Space Flight Center (GSFC) lidar sensitivity level. Remaining line is flight alti- 


tude; arrows represent GSFC lidar pointing direction. 

1 1.5 km and because of cirrus above this region. The two nadir 
profiles below 3 km are similar, particularly across the strong 
cloud layer at about 0.9 km. 

The intercomparison results below 6.5 km for the two wave- 
length regions show remarkably good overall agreement, 
within 5-10%, where aerosol temporal variability was shown 


to be minimal. The precision of this comparison is far better 
than was possible during the transit flights, where modeled P 
at flight level could only be compared with measured P within 
1-2 km altitude from the aircraft [Cutten el a 1996]. In fact, 
the only anomalies in the comparisons discussed so far for P l 06 
data are the constant discrepancy in the FSSP P l06 data at a 
given altitude and the change in magnitude of the discrepancy 



Figure 4. Backscatter profiles near Japan based on 1.06-jim GSFC profiles for several time periods and 
beam orientations (solid, dashed, and dash-dotted lines), GSFC data averaged over ±0.5 km of the flight al- 
titudes listed (vertical bars), and aerosol backscatter modeled at flight level from FSSP (circles) and LOPC 
(diamonds) sensors. Horizontal bars and associated symbols show range of data and averages for each sen- 
sor. Bars at each altitude have been slightly displaced vertically for clarity. 
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Figure 5. Time series of measured and modeled 9.25-pm aerosol backscatter coefficient for selected periods 
from Japan flight on June 3, 1990. Jet Propulsion Laboratory (JPL) data are averaged over ±0.5 km of two 
flight altitudes studied. Plot includes flight level aerosol backscatter (P) data modeled from the FSSP 
(triangles) and LOPC (squares) sensors, and measured P using 9.1 -pm CW Doppler lidar (dotted and dashed 
lines). Solid and dashed horizontal lines at lxl0' n and 6x10 12 m~ 5 sr 1 represent sensitivity levels for JPL 
(near aircraft) and CW lidars, respectively. Remaining solid line is flight altitude; arrows represent JPL Il- 
dar pointing directions. 


with altitude. In cleaner regions of the troposphere (where P 106 
< I0' 8 m' 1 sr J ), smaller particles (<0.5 pm) contribute most to 
p [ Clarke , 1993]. Since the FSSP is unable to sample these 
small particles adequately, this will lead to a considerable un- 
derestimate of P, (}A , providing an explanation for the FSSP dis- 
crepancy seen at higher altitudes. Both show variability in 
modeled P OPCs in the 9.25-pm region, particularly at 9.8-km 
level, which may be the result of small sampling volumes. 


However, the mean FSSP-modeled and CW lidar P (the only li- 
dar data available at the higher altitudes) are similar. This may 
be due to an absence of a significant number of large particles 
(mentioned above), resulting in most of the infrared P origi- 
nating from particles less than 1-pm diameter. Under these 
conditions, modeled phased on FSSP counting statistics will 
not be greatly affected by the problem of the FSSP undersam- 
pling the large particles [Cutten et al . , 1996], but instead will 



Figure 6. Same as Figure 4, but for 9.25-pm wavelength; includes P measured by 9.1 1 -pm CW Doppler li- 
dar (squares). Note no JPL profile for 0242-0255 UTC. 
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be dominaied by counting statistics in the 0.35- to L0-pxn di- 
ameter range. 

4.2. Hawaii Flight 

For the second P intercomparison from remote and in situ 
sensors, complete sets of lidar P data were used for the lower 
two levels (6.5 and 3.2 km) around the big island of Hawaii. At 
the 9.4- and I2.5-km flight levels the GSFC lidar was pointed 
in the zenith direction, while the JPL lidar data contained nu- 
merous dropouts owing to weak backscatter. 

4.2.1. 1.06 pm. Figure 7 shows p l06 data as a time series 
plot for two altitudes. The FSSP and GSFC lidar at 4-km alti- 
tude compare well, with the spread of data from each sensor 
confined to less than an order of magnitude. Here the p was 
dominated by particles above 1-pm diameter where the FSSP 
had moderate counts. On the other hand, LOPC-modeled P 106 
are about 50% higher. 

At the 6.5-km altitude the FSSP-modeled P l()6 show large 
scatter. An examination of the false color time series plot of 
the FSSP size distribution (not shown) indicated that the parti- 
cle numbers within the size range 1.2- to 3-pm diameter fluctu- 
ated over this period. Particle numbers in this size range 
showed two strong peaks (at 2324 and 2336 UTC) where 
counts were observed in all size bins. The suggestion is that a 
well-defined aerosol layer was encountered with orientation N- 
S, since these times correspond to transects of the same longi- 
tude. (Winds at this altitude were from the north at -15 m s' 1 .) 
The GSFC data for 2350-0000 UTC are a factor of -2.5 above 
the FSSP-modeled P U)6 for the leg at 2300-2312 UTC owing to 
(1) the lower cutoff limit of the FSSP sensor and (2) the GSFC 
lidar sample volume probably detecting more P from particles 
larger than !-|im diameter because of spatial inhomogeneities 
associated with moderately low p. Because of its long sam- 
pling times, the LOPC sensor tended to average out these in- 
homogeneities. However, the GSFC data all lie above the mean 
FSSP-modeled P, ()6 , which may be due to the choice of aerosol 


components since the same composition was assumed as for 
the 4-km level. Figure 8 shows two GSFC lidar nadir profiles 
measured from 6.5- and 9.5-km altitude levels at 2300-2345 
and 2350-0000 UTC, respectively, and averaged over the lime 
periods shown in the figure. Again the agreement between the 
two profiles is very good. This p l06 profile comparison, which 
occurred early in the GLOBE II mission and the Japan profile 
comparison from the latter part of the mission, indicates the 
GSFC lidar was performing consistently during the GLOBE II 
experiment. 

4.2.2. 9.25 pm. Figure 9 shows a time series plot of JPL li- 
dar P data and P from each of the flight level sensors at two al- 
titudes (4 and 6.5 km). Agreement to within a factor of -3 or 
less exists at 4 km among the CW lidar and FSSP and LOPC 
sensors. However, the JPL lidar p data were nearly an order of 
magnitude lower around 2340 UTC and over the period 0035- 
0120 UTC. At 6.5 km the CW lidar and FSSP and LOPC sensors 
all show large fluctuations in p. This arose from transecting a 
N-S aerosol plume at 154°W on two sides of the closed flight 
path (see the 1.06-pm discussion above). The FSSP sensor and 
CW lidar fluctuations are to some degree correlated, although 
the amplitudes of the CW lidar fluctuations are smaller by over 
an order of magnitude. The LOPC sensor exhibits smaller fluc- 
tuations but at higher levels. The JPL P measured 1.5 hours 
later at this altitude exhibits larger fluctuations compared with 
4.0 km but much smaller compared with the FSSP. 

The P intercomparison between the lidar and OPC sensors 
in this flight did not show as good an agreement for the flight 
near Japan. More p fluctuations were evident partly due to 
aerosol temporal variations at a given altitude, particularly at 
6.5 km. However, the best agreement seen at 4 km (between 
CW lidar, and both OPCs) supports the previous conclusion in 
section 4.1.1 that the best agreement between lidars and OPCs 
is for moderate P values (>10' 9 m I sr ') 

To summarize, comparisons of measured and modeled P dur- 
ing the two flights showed differing results for the two wave- 



Time (UTC) 

Figure 7. Same as Figure 3 for Hawaii flight on May 18-19, 1990. 
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Pi. 06 ^ 9.25 


Figure 8. Two 1.06- and one 9.25-pm aerosol backscatter profiles averaged over time periods during Ha- 
waii flight and p l t)( / P 925 profile averaged over 2305-2345 UTC. 


lengths. The best agreement at 1.06-pm wavelength occurred the GSFC lidar sampling volume, which in the presence of 
when p i()6 values were >10* m 1 sr' 1 . However, the smaller aerosol inhomogeneity can lead to incomplete sampling and 

FSSP- modeled (J 1<U6 in this p range is partially the result of the consequently more scatter in modeled P u ». The calibrated 

lower diameter cutoff (0.35 Jim), which will cause p l06 to be GSFC lidar, which measured p from particles over a large size 

underestimated by as much as 20%. In tropospheric regions range, is less likely to be affected by small spatial inho- 

where P j0A <10 K m'* sr'\ it appears the P was dominated by mogeneities and will likely show less P l{)6 variation. Agrce- 

particles near or below the FSSP lower cutoff diameter. Fur- ment between measured and modeled P at 9.1 1- and 9.25-jiin 

thermore, both OPC sampling volumes arc much smaller than wavelengths was best for P values >4x10 10 m 1 sr 1 which gen- 



Time (UTC) 

Figure 9. Same as Figure 5 but for Hawaii flight on May 18-19, 1990. 




CUTTEN ET AL: TNTERCOMP ARISON OF PULSED LIDAR DATA 


19,659 


era] Iy occurred below 7-km altitude at both measurement loca- 
tions. For p values <10 10 m 1 sr' 1 , agreement was more diffi- 
cult to ascertain because of higher aerosol variability. 

5. Aerosol Composition 

This section will examine the aerosol composition prevail- 
ing at the two locations using measured p, p ratio, LOPC- 
modelcd p, and CN data. A list of p ratios for several different 
aerosol compositions and wavelengths is given in Table 2, 
taken from theoretical calculations based on single mode log- 
normal size distributions [S rivastava et al 1992]. Here 
P‘>.i/Pio..w provides data on particle composition [Srivastava 
et al 1995]. In low P regimes, P 9 M /p lusy near 10 is due to an 
ammoniated aerosol comprised of ammonium sulfate 
[Srivastava etcil., 1997], while p 9 M /P, osy for partially am- 
moniated aerosol (consisting oi a mixture of ammonium sul- 
fate, sulfuric acid, and dust) has a range of 2-10. However, us- 
ing infrared wavelengths to compute infrared P ratios poses a 
problem because large refractive index variations can occur 
over these wavelengths f Kent et al., 1983], leading to results 
which can be ambiguous. Note that since (3 0 5 ,/p, {)h are small 
these data do not provide information on chemical composi- 
tion. 

5.1. Japan 

The dominant vertical aerosol feature (section 4.1) present 
during the flight near Japan will be examined for aerosol com- 
position using CN time series data, LOPC-modeled p, 9.11- 
and 10.59-jam p, and corresponding backscalter ratio for the 
whole flight (Figure 10). For flight levels above 8.5-km the 
infrared P was generally below 10* 10 m 1 sr' 1 . A low CN refrac- 
tory component (Figure 10c) implied that sulfuric acid domi- 
nated the small-particle regime (<0. 1 -pm diameter). Best 
agreement between the LOPC-modeled and lidar-measured p 
data was found using sulfuric acid composition to compute the 
LOPC-modeled p. The few p 9 n /p, 0V) data above 8.5 km 
(Figure I0d) indicate a ratio around 2, which is also consistent 
with sulfuric acid particles in a low P regime [Srivastava et al., 
1995]. A moderate agreement is noted between the CW lidar 
and FSSP- modeled P at the 11.5-km level, which is also indi- 
cated by the good agreement between the average value data 
(Figure 6). Since P at this altitude was due to moderate concen- 
trations ot sulfuric acid particles with upper diameter limit bc- 


Table 2. Calculated Backscatter Coefficient Ratios for Several 
Wavelength Pairs, Aerosol Compositions, and Geometric 
Mean Radius (r.) Values 



V 

pm 

Sea Salt 

Dust 

Ammonium 

Sulfate 

Sulfuric 

Acid 

P 0 . 53 /P rots 

0.15 

* 

2.5 

3.5 

2.5 


0.6 

2 

2 

2.5 

2 

Pi (*/ P‘J 25 

0.15 

- 

225 

20 

65 


0.6 

25 

5 

2.5 

5 

P'J l/PlO.59 

0.15 

- 

2 

13 

2.5 


0.6 

2-3 

2 

10 

1.5 


N (> - 1 particle cm' 3 and = 1.5 were assumed. 


low 1 jim, the FSSP count statistics are adequate to estimate P 
from these size distributions. At flight levels below 8 km the 
aerosol composition contained a much larger refractory CN 
component (Figure 10c). Both the measured and modeled P 
were greater than ~10‘ 8 m 1 sr' 1 , with P v .,,/p l0 . w remaining less 
than two. This indicates a composition dominated by sea salt 
and/or dust-type aerosol [Srivastava et al ., 1995]. However, 
Plate 1 indicates a plume of aerosol extending upward from the 
MBL to -8.5 km that is expected to be composed predomi- 
nantly of sea salt. FSSP- and LOPC-modeled and CW lidar- 
measured P all show a dip at the beginning of the 6.5-km 
flight level. This indicates that a region was encountered with 
more volatile particles, most likely composed of sulfuric acid 
since p 9 n /p i(JV , remains below two. The large drop in nonre- 
fractory CN at 3.5 km just after 0500 UTC appears to be corre- 
lated with an order of magnitude change in p as measured by 
the CW lidar. This indicates ‘a change in composition, where 
sea salt or other refractory material may have been partially 
replaced by a volatile aerosol. The CW P and the JPL colorized 
P image (not shown) at 3.5 km showed a cleaner region where 
P in each case was an order of magnitude lower. Since the total 
CN concentration increased, the change must be due to an in- 
crease in small sulfuric acid particles, as p (J ,,/p,,, ^ remained 
below two. It is noted that the FSSP-modcled P (near 0500 
UTC in Figure 10b) did not respond nearly as strongly to these 
changes. At 8.5-km, P returned to values below 10' ,n m 1 sr' 1 , 
and p«, n /P MiSlJ exhibited a greater range of fluctuation 
«2 - -3.5). 

An interesting feature was noted in l he p plots shown in 
Figure 10b, over the periods 0315-0400 UTC and 0615-0700 
UTC at 9.8 and 8.5 km, respectively. During these periods the 
300-s CW lidar p, UJ showed only small fluctuations about a 
mean value. Likewise, the JPL lidar p, J35 , averaged 1 to 2 km 
from the aircraft, exhibited a similar behavior. However, the 
300-s FSSP- and much longer averaged LOPC-modeled p data 
show large variations over these two periods. It is noted that 
since the FSSP-modeled p is below 10 10 m 1 sr' 1 , it could be up 
to a factor of 3 loo low [Cullen et al., 1996] than the measured 
P, because of inadequate count statistics. This feature was also 
observed al the 12.6-km flight level in the Hawaii flight and in 
several of the other GLOBE transit flights, all in the infrared p 
range of 10 10 - 10 11 m' 1 sr' 1 . Furthermore, a significant varia- 
tion occurred in the refractory CN over the 0615-0700 UTC pe- 
riod which is correlated with the FSSP- modeled p variation. 

The P 9 i,/p, os «, plot (Figure lOd) exhibits variation over this 
period but is poorly correlated with FSSP-modeled p. The near- 
steady CW lidar p at both wavelengths is a result of their much 
larger sampling volumes. The lidars arc sampling many more 
large particles (>I jam) which contribute to a significant frac- 
tion of the P compared with (he much smaller sample volume 
of the OPCs. If the composition of these large particles was 
constant over (he period, then little variation would result in 
the lidar p. On the other hand, variations seen in the OPC data 
arc probably due to sampling aerosols <l-jam diameter whose 
composition varied along the flight path. For example, during 
the 0615-0700 UTC period, a small sulfuric acid aerosol com- 
ponent was initially encountered, which changed to a sulfuric 
acid-dust mixture that left a nonvolatile component after heal- 
ing and which finally reverted to mostly sulfuric acid. 

The analysis for the Japan flight has revealed two interest- 
ing features. First, the aerosol plume, evident in the 1.06-jim 
colorized p plot (Plate I), contained predominantly MBL sea 
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Figure 10. Flight level measurements of (a) condensation nu- 
clei concentration, (b) measured and modeled 9.11-p.m aerosol 
backscatter (horizontal line represents CW lidar sensitivity 
level), (c) refractory condensation nuclei concentration, and 
(d) ratio ofp y .„/P, 0.59 aerosol backscatter, for entire flight near 
Japan. 


salt aerosol which had risen to -8.5 km. The lower CN concen- 
trations (Figure 10a) at 3.2- and 6.5-km flight levels below 8 
km indicate the presence of an aerosol predominantly from the 
MBL. The aerosol plume may have also contained sulfuric acid 
particles (<i-|im diameter) that originated from pollution 
sources in Japan. One layer observed at 3.5 km was probably 
the result of concentrations ol pollution-type particles. Above 
the aerosol plume (-9 km) all measured and modeled (3 data, as 
well as P ratio data, indicate a free tropospheric sulfuric acid 
aerosol which is normally found at these altitudes. Second, it 
was found that in certain aerosol regimes, both of the in-flight 
OPC sensors showed a widely varying, modeled p„ 25 , while the 
lidar-measured p 925 showed minimal variation. This can be ex- 
plained if the large particles (>2-pm diameter), which the lidars 
would sample adequately, did not change composition and the 
smaller particles did. Thus a coarse form of composition seg- 
regation occurred when using sensors with very small and very 
large sampling volumes to inter P values. This augments the 
point that a fully representative infrared P whose real value is 
<10' ,n m ' sr' 1 can only be obtained if the sensor with a large 
sampling volume is used [Cutten et ai , 1996]. 

5.2, Hawaii 

The microphysics analysis was restricted because of less 
comprehensive coverage of the lidar data compared with the 
Japan flight. However, OPC-derived p data combined with one 
ratio (pj 06 /p y 2S ) profile (Figure 8) provide some insight into 
the aerosol microphysics prevailing at the time in the lower 
altitudes. It should be noted the two averaged P, £ * profiles in 
Figure 8 (one over 40 min. the other over 10 min) compared 
very well below 5 km, indicating the longer averaging peiiod 
encompassing all three legs did not produce significant differ- 
ences from the first leg of the flight. The nadir p, U6 /P<>.25 P ro ~ 
file shows values ranging from 25 to about 110 where the 
lower p, ()6 /P«j .2 5 values correspond to a wavelength depend- 
ence exponent of ~1.5, implying the presence of large dust 
and/or ammoniated particles over the altitude range 3.6-4. 1 
km. In addition, the CN data measured at 4-km flight level be- 
tween 2200 and 2300 UTC indicate that a strong refractory 
component was present for this size range. In Figure 9 LOPC- 
modeled p,, n showed best agreement with CW lidar p„ , t using a 
composite aerosol mixture based on modeling by Srivastavci et 
ai. [1997] which was modified with a hygroscopic refractory 
component representative of, for example, sulfuric acid coated 
dust. LOPC- derived P UI6 were slightly higher than the FSSP- 
modeled P, <16 for the same modeled composition. It is noted 
that when the hygroscopic dust component was replaced with a 
nonhygroscopic dust, the calculations yield a noticeably lower 
LOPC-modeled P 0 At 2-km altitude the P ratio curve (Figure 
8) indicates a decrease in the concentration of large aerosol 
particles with height and a possible change in composition to 
a partially ammoniated aerosol; there is only a small decrease 
in P 106 . It should be noted that spatial resolution of the data is 
not a factor, since vertical resolution ot the GSFC lidar is half 
that of JPL lidar. The wavelength exponent for aerosol P in 
this layer is around 2.2. Below 2 km the ratio fluctuates around 
a value of about 55 (wavelength exponent of -1.95), where 
both P l( *and p, J2 < are high. In this case it would appear that a 
refractory material again was responsible for the high p. Thus, 
from near the surface to 5 km, there is a change in large parti- 
cle numbers and/or particle composition; this behavior may 
have been caused by aerosol plume emanating from volcanic 
activity that was occurring at that lime on the big island of 
Hawaii. 
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6. Conclusions 

The main conclusion from the intercomparison is that all 
measured and modeled p data showed reasonable agreement at 
moderate P levels but lesser agreement in low P conditions 
which often occur in the middle and upper troposphere. In the 
latter case the pulsed lidar data showed the largest disagreement 
with the flight level sensor data, whereas differences within 
the in-flight sensor data were generally smaller. This signifi- 
cant difference in P between the OPCs and pulsed lidars may be 
due to the large differences in sample volumes, the sampling 
volume factor being as much as I0 8 between the pulsed lidars 
and the FSSP sensor. Thus the lidar P will be more representa- 
tive for both wavelength regions because (1) the lidars sample 
the larger particles (which contribute most to p for the 9.25- 
pm region) more frequently, and (2) there is no comparable 
cutoff in the submicrometer particle size range (which the 
OPCs exhibit) that affect the 1.06-pm region. Limited data 
from the upper tropospheric region did not indicate that the 
FSSP-modeled infrared P was underestimated when 
P <10 ,f} m‘* sr 1 \Cutten et al , 1996]. 

This study has also revealed that a suite of airborne instru- 
ments, like the one described here, can provide detailed infor- 
mation on the prevailing aerosol. The two distinct aerosol re- 
gimes identified in the Japan flight, along with smaller 
changes within these regimes, attest to this. Using lidar p ra- 
tio data for several wavelength pairs can reduce the ambiguity 
in determining aerosol size ranges responsible for P at a par- 
ticular wavelength. Furthermore, P ratios, when combined 
with theoretical calculations using aerosol constituents com- 
monly found in the tropospheric aerosol, can allow further in- 
formation to be retrieved remotely about the composition 
[Srivastava et al . , 1997]. 

Finally, this study has highlighted the usefulness of corre- 
lating data from measured P ratios with CN refractory and non- 
refractory particle numbers to derive the aerosol composition. 
It also provided evidence that aerosol chemical composition 
can be segregated based on sample volume size. Large lidar 
sample volumes can sample large particles much more effi- 
ciently; when the particle composition differs from that of 
smaller particles, the lidar measurements can give a different P 
from that modeled using the OPC size distributions. 
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